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Abstract Long-exposure spectro-polarimetry in the near-infrared is a preferred 
method to measure the magnetic field and other physical properties of solar 
prominences. In the past, it has been very difficult to observe prominences in this 
way with sufficient spatial resolution to fully understand their dynamical proper¬ 
ties. Solar prominences contain highly transient structures, visible only at small 
spatial scales; hence they must be observed at sub-arcsecond resolution, with 
a high temporal cadence. An adaptive optics (AO) system capable of directly 
locking-on to prominence structure away from the solar limb has the potential 
to allow for diffraction-limited spectro-polarimetry of solar prominences. In this 
paper, the performance of the off-limb AO system and its expected performance, 
at the desired science wavelength Can 8542 A, are shown. 

Keywords: Adaptive optics, Solar prominences 


1. Introduction 


Solar prominences consist of relatively cool, dense plasma, which is suspended 
above the solar surface. Filaments and prominences represent the same phe¬ 
nomenon, viewed differently ( |Tandberg-Hanssen, 1995 Labrosse et ai, 2010| [MacKa y et al ., 2010| ).[ 
The prominence itself forms within or directly beneath a magnetic flux rope 
dMacKay et al., 2010||Berger, 2014[ ). A flux rope is a long, twisted magnetic field 
region, which is nearly force-free (Tandberg-Hanssen, 1995; Ma cKay et al., 2010[ ).[ 

Such a flux rope may form either from the shearing of magnetic arcades, or it may 
emerge from below the photosphere, already formed dMacKay et ai, 2010p . One 
main reason that the study of solar prominence should be considered important is 
that they are often involved in coronal mass ejections (CMEs). CMEs occur when 
a flux rope becomes unstable, e.g. when there is emerging magnetic flux, coining 
from below (Chen, 2011). Since flux ropes very often contain solar prominences, 
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prominences may show instabilities that could be used to predict an imminent 
CME ( |Chen, 201 ip . 

Despite continuing study of solar prominences, there is still much that is 
unknown about their physical and magnetic structure at small spatial scales 
dLabrossc et al., 2010} |MacKay et al., 2010) |Bcrgcr, 2014[ ). There are ambigu¬ 
ities related to the broadening of spectral lines of prominences that are at¬ 


tributed to unresolved fine structure (Labrosse et ai, 2010). Thus the ability 


to measure spectra of solar prominences at very fine spatial scales is neces¬ 
sary to the understanding of solar prominences. Since solar prominence plasma 
interacts with the magnetic field in which it resides, an understanding of promi¬ 
nence dynamics can only be understood by the inversion of spectro-polarimetric 
data ( |MacKay et al, 2010[ ). There are many ambiguities about prominence be¬ 
havior, with respect to their magnetic fields, at the small spatial scales be¬ 
cause these structures are unresolved dMacKay et al ., 2010 1. Space-based in¬ 
struments are capable of imaging solar prominences at very fine spatial scales 
dBerger et al, 20TTj ). However, to the knowledge of the authors, none has been 
launched that can perform spectro-polarimetry on a solar prominence, which is 
necessary for the understanding of prominence magnetic fields dTandberg-Hanssen, 1995| [ 

|MacKay et al ., 2010| ). 

Ground-based telescopes have instruments which are capable of taking spec¬ 
tral and spectro-polarimetric data, for example the Interferometric Bldimen- 
sional Spectrometer (IBIS: ICavallinil [2006D and the Facility Infrared Spectrome¬ 
ter (FIRS jJaeggli et ~aZT| I2010]l on the Dunn Solar Telescope (DST) in New Mex¬ 
ico. Current solar adaptive optics (AO) systems are designed to utilize broad¬ 
band light and are thus confined to locking onto structures on the disk of the Sun; 
prominences are invisible in broadband light (Rimmele and Marino, 20ll|. This 
limits the usefulness of current solar AO systems. They can only be used to cor¬ 
rect prominence images, and hence provide high resolution spectral or spectro- 
polarimetric data when there is a pore or other dark feature directly adjacent to 

the part of the limb near the prominence, as was done by |Orozco Suarez, Ascnsio Ramos, and Trujillo Bueno|| 

(120131) . These data, however, are fundamentally limited in resolution because AO 

systems can only provide their best correction extremely close to the point upon 

which they are locked ( [Hardy, 1998[ ). Only a purpose-built AO system that can 

directly lock onto solar prominence structure can allow for spectroscopic and 

spectro-polarimetric data at the diffraction limit of the telescope; thus allowing 

for an increase in humankind’s understanding of solar prominences. 

The authors are the first group to construct an off-limb solar AO system 
mainly because of the difficulty in measuring the incoming wavefront using only 
light from solar prominences. The main issue is photon flux. Ha light from solar 
prominences was chosen for wavefront sensor (WFS) of the off-limb solar AO 
system because prominences emit very brightly at this wavelength, relative to 
other spectral lines. Even so, the required 0.5 - 0.7 A filter bandwidth transmits 
very little flux. Thus, it is necessary to use all available Ha light for the WFS. An 
entire AO system was designed and optimized around this WFS. The difficulty 
of measuring the wavefront, coupled with the need to design an entirely new AO 
system, is the reason why the authors are the first to build an off-limb solar AO 
system. (For a review of specific problems of solar AO in general, see the earlier, 
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comprehensive reviews of IRimmele and R adick (1998) and IRimmele and Marinol 

(HsnD-) 

A correlating Shack-Hartmann wavefront sensor (SHWFS) was chosen for 
the off-limb solar AO system, due to its well known characteristics, though it is 
not the only type of WFS that was considered ( |Taylor et al ., 2012| ). A standard 
SHWFS uses point source images directly, i.e. from a star or laser, to measure 
the wavefront. It does this by breaking the beam of light into small pieces, each 
representing a portion of the telescope aperture. The position of each star image 
directly corresponds to the wavefront slope at each position in the telescope 
aperture. A map of the wavefront is created from these slope measurements 
(| Hardy, 1998 |); see Figure [5] Correlating SHWFS have been used in solar AO 
for several years fRimmcl c and Radick, 1998] |Rimmele and Marino, 20TT| ). The 
difference between a correlating SHWFS and a conventional one is that each 
sub-aperture image is cross-correlated with a reference image, typically one sub¬ 
aperture image from a single frame. The maximum of each cross-correlation func¬ 
tion shifts with the image motion, and thus the changing wavefront slope in each 
of the SHWFS sub-apertures ( Rimmele and Radick, 1998l|Rimmele and Marino, 2011[ ).| 

In this paper the results of some of the early tests of a closed-loop AO system 
on the DST are shown, based on a correlating SHWFS. In Section[2l |Taylor et al.\ 

(120181) is reviewed, which details the performance which was expected to be 
achieved from the off-limb solar AO system. In Section 01 the setup of the 
experiment which was used to verify the performance predictions is detailed. 

In Section 01 the results of these experiment are enumerated and it is shown 
how they relate to the predictions in Section [2] 


2. Predicted Performance 


As was stated in Section [T| there is very little light available for the SHWFS. 
Indeed, there is so little flux available, even in Ha light, a custom made Ha filter 
with very high transmission was used; see Figure [1] This same type of filter will 
be used on the Visible Broadband Imager on the Daniel K. Inoue Solar Telescope 
(DKIST, formerly ATST; IRimmclc et ail 12014 ). A Dalsa Falcon VGA300 HG 
machine-vision camera was selected, due to its high sensitivity, moderate read- 
noise, high speed, and low cost, coupled with an in-house lenslet array, to make 
the SHWFS. To achieve the required speed, only a small portion of the camera 
chip (160 x 160 pixels) was read-out. 

Using the high transmission filter and sensitive camera, the SHWFS is still 
photon starved; thus it is necessary to make trade-offs in design parameters. If the 
wavefront were finely sampled, by dividing the SHWFS field of view into many 
sub-apertures, a great deal of light would be lost and the SHWFS would no longer 
be able to sense the wavefront at a sufficiently high frame-rate to allow for AO 
correction. Therefore, fewer sub-apertures were used for this SHWFS than for the 


SHWFS used for on-disk observation at the DST (Rim mele an d Ra dick , 1998). 
A large field of view (FOV) per sub-aperture is also needed, to allow for the 
tracking of solar prominence features that are large in comparison to features 
seen on the disk ( [Taylor et al., 20l2| . Finally, the final pixel scale of the sub¬ 
aperture images must be chosen to balance resolution and image brightness 
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Figure 1. Filter transmission profile of the custom Ha filter. 


QMichau, Rousset, and Fontanella, 1993[ ). The number of sub-apertures used, 
their FOV, and the final pixel scale can be determined analytically, as well as 
directly measured ( |Scharmer et al., 2003' ). To verify that an off-limb solar AO 
system is possible and practical, the performance of a particularly promising 
SHWFS configuration was modeled in detail. 

2.1. WFS Noise from Indirect Methods 

It is necessary to measure the noise created by the SHWFS and to estimate 
noise from the various other sources that are present in a finished AO system, 
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in order to predict the ultimate performance of an off-limb solar AO system 
( [Hardy, 1998) . 

A SHWFS was setup with a 5 x 5 array of sub-apertures (see Figure O, each 
with an FOV of 30" and an image scale of 0.80" per pixel. The image scale 
and FOV were measured using a resolution target. This camera was run with a 
frame rate of 900Hz and exposure times of 900ps. The functional form for noise 
generated by a correlating SHWFS is given by IRimmele and Marinol (1201111 . 
However this function was explicitly stated to be defined when the Nyquist 
sampling criterion is satisfied, which this system cannot satisfy, due to low light 
levels. |Michau, Rousset, and Fontanellaj (119931) give a more complete version of 
this formula: 


5m 2 a 2 {dp) 2 

= riR (/A) 2 


[waves 2 ], 


( 1 ) 


where a 2 is the total noise and m is the full-width at half-maximum (FWHM) 
of the cross-correlation (CC) peak, from the sub-aperture images (see Figure 
The quantity er{( is the detector noise, which is assumed to be read-noise 
plus photon noise. The read noise was measured to be approximately 68 elec¬ 
trons per pixel (see |Barry and Burnell] 12000 . for details). The photon noise 
was assumed to be Poisson noise: cr 2 hoton = ^SadU) where g is the camera 
gain and Sadu is the average pixel value, in analog-to-digital units (ADUs) 
dBarry and Burnell, 2000 ). The term n 2 of is the energy content of the image, 
d is the sub-aperture size, p is the physical size of the pixels, / is the effective 
focal length of the optical system (telescope plus bench optics), and A is the 
wavelength at which the measurements are made, 6563 A. Note that 


{dp) 2 {d0) 2 


( 2 ) 


where 9 is the image scale, in radians per pixel, and d is 13.07cm (see Figure [2]). 

To measure m, the CC peaks for the 5x5 array were used (Rim mele and Marino, 20 11]).[ 
Following the method of Mich au, Rousse t, and Fontanella (]1993j) , a two-dimensional] 

(2D) Gaussian function was fit to the central portion of each CC peak. The 
value of m is the average FWHM of all the CC peaks.The detector noise is again 
denoted as fib. 

The energy of the image is 


n 


2 2 
„ ( T: 


Ei 5 ^i 2 ’ 


( 3 ) 


where Sj and S are the signal per pixel and the average pixel value, measured 
in electrons ( |Michau, Rousset, and Fontanella, 1993[ ). 

The result of Equation HD was calculated, for each of the 2000 frames, at each 
exposure time, and then averaged over those 2000 frames. To make the smooth 
curves, as shown in FigureEJ one should note that |Michau, Rousset, and Fontanella] 
(119931) used an approximation for an image of a given contrast, using the above 
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Figure 2. The pupil layout for the 5x5 grid. The outer circle is the telescope pupil, 76.2cm 
in diameter. 



Figure 3. SHWFS images, using the Dalsa Falcon VGA300 HG camera. On the left are 
prominence images which have been dark-subtracted and flat-fielded. On the right are the 
cross-correlation peaks. 


notation: 


5 m 2 cr 2 ( d9 ) 2 

n^CgS adu 


[waves 2 ], 


( 4 ) 


where C is a factor related to the image contrast and n 2 is the pixel area of 
each sub-aperture image. The value of C is found by dividing Equation © by 
Equation ©. Thus 


C = 


5<Sadu 


( 5 ) 


for an exposure time of 900ps. The results were then plotted by scaling a\ and 
u Photon with exposure time. 
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Averaged Tilt Power Spectrum 



Averaged Tilt Power Spectrum 



Figure 4. Noise power spectrum for data taken at 900 Hz, for moderate (top panel) and high 
(bottom panel) contrast prominences. 


2.2. WFS Error from Telemetry 


In order to verify the above error measurements, the error in the SHWFS was 
found by direct measurement from SHWFS telemetry. This was done by tak¬ 
ing the power spectrum of the total pixel shifts as measured, that is, shift = 
\J £g hi f t + 2/shift ■ The power spectrum should be taken for data with a zero mean, 
as stated by IMarinol (120071) . So the average shift was subtracted from each 
member of the shift array, for each data set. 

Using the method of IMarinol (12007|) . assuming white noise, we obtain 
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( 6 ) 


N 

ct wfs = E Noise level. 

1 

Hhere the ’Noise level’ is the point where the power spectrum becomes flat and 
N is the total number of exposures in each time series (see Figure a. This 
becomes 


°wfs = Noise Level • N, (7) 

which yields the noise variance in terms of pixels 2 . To find the RMS noise in 
terms of nanometers, the following is applied, 


^WFS nm = ct W fs^/(3600-57.3). (8) 

Here 6 is the pixel scale, 0.82” per pixel, and d is the sub-aperture diameter, 
13.07 cm. 

To find the error in radians, we use 


0"WFS rad = 27T(T W FS nm/A, (9) 

with A being 6563 A. The ratio between ctwfs rad and a x , as found in Section ETT1 
was calculated for each exposure time where ctwfs rad was defined. The value of 
a x was then scaled by the average of these ratios, for each prominence. Equation 
(|4|) was then fit to the scaled a Xl as shown in Figure [5] (The error due to the 
rate at which a SHWFS can measure the wavefront is also shown, as explained 
in Section 12.31 1 Comparing the result from Equation dTJ) with the telemetry 
measurement, Equation m is found to overestimate the noise by about a factor 
of 4. 

Figure [5] shows the importance of balancing the speed of the wavefront mea¬ 
surements; the shorter the exposures, and thus the faster the wavefront mea¬ 
surements, the lower the noise from the atmosphere, but the higher the noise 
from the SHWFS itself. The way in which this impacts the quality of the final 
image in will be shown in Section 12.31 

2.3. Strehl Ratio 


The quality of an image is often measured in terms of its Strehl ratio. The Strehl 
ratio is the ratio of the intensity of the peak of a point source image measured 
by an imperfect optical system against the intensity of that same object as 
measured by a perfect optical system. A perfect optical system will have a Strehl 
ratio of 1. Any imperfection in the optics or caused by atmospheric distortion 
will lower the Strehl ratio (Hardy, 1998 Born and Wol f, 1999 ). A Strehl ratio 


above 0.8 is considered diffraction-limited ( |Born and Wolf, 1999| ). However, an 
image still contains all of the information of a diffraction-limited one, only 
at a lower signal-to-noise ratio, until the Strehl ratio drops below about 0.1 
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QHardy, 1998||Lukin and Fortes, 1998| . However, the higher the achieved Strehl 
ratio, the higher the signal-to-noise ratio of the data. 

In order to determine the Strehl ratio that one can expect to achieve from an 
AO system using this SHWFS, it is necessary to account for the other primary 
sources of error: time delay, fitting and aliasing errors. For all of the following 
error estimates, a Fried parameter ( |Fried, 1965P , vq, of 10cm, at 5000 A was 
assumed (See |Hardy|(11998jl for a more detailed discussion of AO error sources.). 

To determine the error due to time delay, the functional form quoted by [Hardy] 

(|1998l) was used: 

4d = 28.4(rJ g ) 5 / 3 , (10) 

where r s is the time delay, and f g is the Greenwood frequency (/ g « 0.427//- for 
a single layer atmosphere model, which was assumed for this study), v being the 
wind speed ( |Hardy, 1998[ ). This error only takes the speed of the camera into 
account. 

To find r s , the delay was assumed to be equal to the exposure time plus a 
constant, pessimistic 500ps delay, for wavefront reconstruction plus the camera 
readout, etc. This yields r s = 500ps + exposure time. This is similar to the 
method used by |Hardy| (119981) . Three curves were calculated for ro of 10cm 
at 5000 A, which is 13.8cm at 6563 A, using the formula r' 0 = ro (A/5000) 6 / 5 
dHardy, 1998p . Wind speeds of 5, 10, and 15m s -1 were chosen. Thus, the 
Greenwood frequencies were 15Hz, 31Hz, and 46Hz (see Figure [5]). The formula 
for r s above does not take into account the driving frequency of the deformable 
mirror (DM), which was unknown when this analysis was performed. 

The forms of fitting and aliasing errors were used, as described bv IRimmele and Marinoji 
(120111) . where the fitting error is 


/d\ 5/3 

<7p = 0.28 f — J ■ (11) 

The aliasing error is 

/ d \ 5/3 

°a = 0-08 y—j , (12) 

where d is the sub-aperture size. 

In Figure [6l the resultant Strehl ratios, given the above errors, are shown. 
The approximation Strehl ratio = e~ was used to find these Strehl ratios 
dHardy, 1998| ). All of the Strehl ratios were calculated at 8542 A. To convert the 
SHWFS and time delay errors at 6563 A, which were measured in radians, the 
authors took note that 1 wave = 27 t radians at a given wavelength, so 


08542 in radians 


6563 

8542 


06563 in radians. 


(13) 


Figure [6] shows that one can expect to see Strehl ratios of between 0.6 and 0.7, 
at 8542 A, given a Fried parameter of 10cm and 10m s 1 wind. 

There are a few important facts that can be gleaned from Figure [6] (1) The 
faster the SHWFS is able to measure the wavefront, the better the achieved 
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Strehl ratio. (2) Using faint, low contrast prominences for sensing the wavefront 
lowers the maximum achievable Strehl ratio. (3) For a fainter prominence, the 
maximum achievable Strehl ratio occurs at a slightly slower frame-rate than for 
a brighter prominence. (4) This SHWFS is limited by hardware to a speed of 
around 900Hz, which corresponds to an exposure time of about 900ps. Given a 
faster camera, slightly better Strehl ratios could be achieved, but the achievable 
Strehl ratio decreases slowly, to the right of its maximum point. 


3. Experimental Setup and Procedure 


The off-limb AO SHWFS, as tested in Section^ was integrated into AO bench of 
the DST. A 97-actuator, Xinetics deformable mirror (DM) and a tilt-tip mirror 
(TTM) were utilized. Wavefront reconstruction is accomplished by utilizing a 
customized implementation of the Kiepenheuer-Institute Adaptive Optics Sys¬ 
tem (KAOS), which was originally coded for the German vacuum tower telescope 
(VTT) and was rewritten to operate the GREGOR telescope ( |Bcrkefeld et ai, 2012| .| 


It was customized to utilize the layout of the DM as well as the layout, exposure 
time, etc., of the SHWFS. This implementation of KAOS is run on an off-the- 
shelf computer, utilizing an Intel i7, quad-core processor. This is interfaced to the 
DM using a Xinetics, Gen III chassis, driving 144 channels. Of the 144 channels, 
97 channels are used for the DM, 2 for the TTM, and the remaining channels 
are unused. 

Tests were performed on Port 2 of the DST. The setup consists of a beam¬ 
splitter, which diverts 95% of the incoming Ha light to the SHWFS, leaving 
the remainder to feed the imaging camera. The rest is fed to the SHWFS. The 
bench optics shown are required to create a 5 x 5 SHWFS array with sub-aperture 
FOV of 30" and pixel scale of 0.82" per pixel (see Figure [7]). For these tests, the 
camera was run at a frame rate of 860Hz. This is slightly slower than the model 
frame rate, but it allowed for a border of a few pixels around the sub-aperture 
images (174 x 174 pixels are now being read-out). This allows the KAOS system 
to compensate for slight misalignments in the optical system. Since the system 
is running at a slightly slower frame rate than used in Section [2j an exposure 
time of 1ms could be used. This lowered the noise coming from the SHWFS, at 
the expense of slightly more time-delay error. 


3.1. Procedure 


The KAOS system is capable of recording data which contains raw measurements 
from the SHWFS as well as the calculated errors in the reconstructed wavefront. 
KAOS records the commands sent to the TTM and the DM. It calculates the 
total atmospheric error in each mode, given the residual error and the commands 
sent to the DM. These data and their caveats are explained below. 

When KAOS reconstructs the wavefront, it does so by utilizing Karhunen- 
Loeve (K-L) modes ( [Noll, 1976] |Dai, 1995| . The degree of correction that can 
be achieved depends upon the number of modes that can be utilized, which in 
turn depends upon the resolution to which the SHWFS can sample the wavefront 
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flNoll, 1976||Hardy, 1998[ ). Even if those K-L modes can be perfectly determined, 
there is a residual error in the wavefront, corresponding to the infinite number 
of K-L modes which were not sensed ( |Noll, 1976] |Dai, 1995] ) . K-L modes are 
ranked in various orders. Each mode in each consecutive order contributes less 
to the total error in the wavefront. Thus, the correction of a few modes can 
correct most of the distortion in the incoming wavefront (Noll, 1976 |Dai, 1995] ) . 
This setup can sense 20 K-L modes. fNollj (U976[) has tabulated the residual error 
expected after correcting 20 Zernike modes, which are similar to K-L modes, but 
with slightly higher residuals errors ( [Noll, 1976} |Dai, 1995] ). 

The other major caveat is the way in which KAOS arrives at its measurements 
of the total wavefront error in each K-L mode. It is done by noting the residual 
error in each mode and determining how much correction is being applied in 
each mode via the TTM and the DM and combining the two. If the pixel 
scale of the SHWFS was measured incorrectly, both the residual and total error 
measurements will be wrong. However, these measurements of the pixel scale 
have been ensured to be accurate. The DM used was a spare at the DST, and 
its neutral state was not exactly flat. This coupled with any static errors in the 
system, which would be corrected by the DM, affects the total wavefront error 
measurement s. 

To compare the measured performance of the off-limb solar AO system with 
the expected performance, shown in Figure [Gj there are a few parameters which 
are of interest: ro, B c \ (the closed-loop bandwidth), and the Strehl ratio of 
the system. If the AO loop is treated like an RC filter, where higher update 
frequencies are attenuated, then B c \ can be defined as the Odb bandwidth. This 
is the highest frequency at which accurate wavefront corrections can be made 
( [Roddicr, 1999[ ). See the explanation of Strehl ratio in Section 12.31 

The value of r' 0 can be calculated from the formula given by : Nolll (11976[) , 


4? = C • (D/r') 5/3 - (14) 

Here D is the diameter of the telescope and <fi 2 is the variance wavefront error, 
in radians, averaged with time. C is a constant which is defined by the expected 
relative variance in a given mode, for a given turbulence model ( [Noll, 1976} 
|Dai, 1995] ). The average variance of the first three K-L modes, after the two 
modes corresponding to tip-tilt, were used, since these have the same approx¬ 
imate variance of 0.023927(-D/rg) 5 / 3 ( [Dai, 1995] ). Also, tip-tilt mode measure¬ 
ments are sensitive to vibrations in the room, etc. The quantity r' 0 was calculated 
by finding the average wavefront error over a series of 10,000 frames, or approx¬ 
imately 12 s. With r' 0 known, it is trivial to find r g. The Strehl ratio can be 
approximated by( |Hardy, 1998| ) 

e -!x[(Ei 0?)+residual]. 

Here JN cj) 2 is the sum of the residual errors for each mode over the time series, 
calculated by KAOS, and the residual is approximated here as 
0.0208(-D/7 q) 5 / 3 ( |Noll, 1976| . This is the residual for Zernike polynomials, but 
it is only slightly more than that for K-L modes ([Noll, 1976] |Dai, 1995]). 
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Finding B c \ is a bit more involved. It is found by taking the power spectrum 
of the residual error for a given mode with the AO system off, and plot it against 
the power spectrum of the same mode with the AO system on (see Figure 0. 
The point at which they first cross is the Odb B c \. The Welch method was used 
to plot the power spectra (jWelch, 1967). 


4. Results 

The telemetry data taken according to the procedure found in Section 13.11 were 
analyzed. These data were taken with the AO system correcting for 20 K-L 
modes. The Odb bandwidth was found to be approximately 38.8Hz (see Figure 
0. This value depends upon the degree of smoothing of the power spectra, so it 
is only a ” ballpark figure”. 

A sample of the r o values and Strehl ratios, which were calculated using the 
above method, is shown in Table 0 The data represent a temporal average over 
a period of 12 s. One can note the large standard deviation among the ro values. 
This is due to the rapidly changing atmosphere. Since total Strehl ratio depends 
upon the Greenwood frequency, which was not measured with the available data, 
it is no surprise that there is variation in measured Strehl ratios, given similar 
ro values. Even given this fact, the below results measure favorably against the 
above predictions, as shown in Figure [6] It should be noted that even with the 
very large values of ro present in some of the data sets, correction of 20 K-L 
modes provides higher Strehl ratios than tip-tilt correction alone. This can be 
seen by using Equation m, 

4> 2 = C- {D/r' 0 f'\ (16) 

If only tip-tilt errors are corrected, the constant C is equal to 0.134 QDai, 1995| ). 
This yields a residual error at 5000 A, at ro = 17.5cm, of 1.56 radians 2 . This, 
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Table 1. A sampling of ro and Strehl ratios, as calculated in Section [3.11 The standard 
deviations quoted show the spread in the data over 12 s: Seeing is often quite variable. 


rg cm, 5000 A 

Standard deviation 

Strehl, 8542 A 

Standard deviation 

Strehl, 6563 A 

15.9184 

7.5325 

0.7926 

0.0652 

0.6773 

14.9362 

6.9773 

0.7884 

0.0610 

0.6709 

17.4868 

9.2392 

0.8456 

0.0395 

0.7536 

9.1684 

4.7710 

0.6605 

0.0609 

0.4978 

6.4664 

3.1955 

0.4789 

0.0713 

0.2911 

14.0340 

6.7349 

0.7797 

0.0791 

0.6601 

12.1283 

6.4545 

0.7635 

0.0497 

0.6346 

6.3531 

3.2902 
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Figure 7. The experimental setup on Port 2 of the DST. The WFS is at top, with the imaging 
camera at left. An interferometer can be used to verify the flatness of the DM. It is to the right 
of the WFS and is used by diverting all light from the WFS, via a mirror, along the rightmost 
red path. 


Power Spectra 



Figure 8. The power spectra from which B c \ was calculated. The approximate Odb bandwidth 
value was found by interpolation. 
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Figure 9. A graph of Strehl ratio vs. ro value. This takes into account the residual wavefront 
errors plus the expected variance left-over from only correcting a finite number of modes 
(Noll, 1976). Note that for ro ~ 10cm, the Strehl ratio is between 0.6 and 0.7, bettering the 
predictions in Section 2.3 
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Figure 10. The above images were taken on 10 December 2013, between 16:50 and 16:55 
UTC. Each image is the average of 100 frames, with a cadence of 4s —1 . Each series was taken 
at 6563 A. They show what a long exposure would look like during moments of moderate 
seeing, with the AO system totally off (left), only correcting for tip-tilt errors (center), and 
with the AO system correcting for 20 modes, ro varied from between about 15cm and 17cm, 
when these images were taken. (The second line in Table [T1 was taken a few moments after the 
right hand image was taken.) Cross sections of each image are shown by the white line. Pixel 
values for each will be shown in Figure [TT] 

converting to radians 2 at 8542 A, as in Equation ED. and calculating the Strehl 
ratio: 

Strehl ratio = e residual variance , (17) 

yields the maximum possible Strehl ratio at 8542 A, only correcting for tip-tilt 
errors of 0.59. The maximum possible Strehl ratio, given correction of 20 modes 
is 

0.92 using a value of C of 0.0208 ( [Noll, 1976) . This system achieved 0.85. So 
theoretically the improvement is small for excellent seeing, but when ro is 7.7, 
the maximum possible Strehl ratio is 0.12, using only tip-tilt corrections. The 
maximum possible Strehl ratio using 20 K-L inodes is 0.72; this system achieved 
an estimated Strehl ratio of 0.6. 

Also note that when the Strehl ratio is highest, its standard deviation is 
lowest, which points to the stability of the off-limb solar AO system, during 
times of good seeing. The 5000 A ro and 8542 A Strehl ratios from Table [flare 
plotted in Figured] Note that for r 0 around 10cm, the Strehl ratio is between 
0.6 and 0.7, slightly exceeding the above predictions. This may be due to the 
very pessimistic time-delay error estimation. The spread in ro vs. Strehl ratio is 
almost certainly due to the variation in Greenwood frequency, but this is seen 
to be a secondary effect. Recall that in Section HPI it was stated that there are 
error sources that these measurements do not take into account. Therefore, the 
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Image Cross Section, AO Off 




Image Cross Section, AO On 



Pixel Number 

Figure 11. Cross sections taken of each image image in Figure ITOl The increase in sharpness 
for the fully corrected image is depicted by the sharp edges shown in the graph. This indicates 
that much finer details are visible. 


Strehl ratios in Tabic |T] and Figure [9] will almost certainly be higher than the 
true Strehl ratios, as might be measured at the final prominence image. Figure 
1101 depicts a few images which show the dramatic improvement that the off- 
limb solar AO system can provide. The images are each made from a series of 
100 frames with an exposure time of 200ms each. They were taken at a frame 
rate of 4s _1 and span 25s. To better show the resolution improvement made 
in this figure, their cross sections are plotted in Figure [TT] These images and 
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Figure 12. Wavefront error in radians at 8542 A. These data have been smoothed by convolv- 
ing them with a Gaussian kernel, with a standard deviation of 15. The total error calculation 
becomes very noisy, once the TTM is activated, but the total, subtracting tip-tilt errors stays 
well behaved, even when the AO is fully activated. Since the reconstructed tip-tilt errors were 
not used for any calculations, this is of little consequence. 


cross sections show that even during good seeing, long integrations, which will 
be required for spectro-polarimetry, are vastly clearer with the off-limb solar AO 
system fully operating. During very bad seeing, the tip-tilt mirror can be used 
alone, which is still an improvement over limb tracking, the standard method 
for tracking prominences on the DST (Kevin Reardon, private communication). 

In Figure [T21 the improvement in image quality given by the off-limb solar AO 
system is shown. For the first 3000 frames, the AO is applying no correction, 
the residual error is the same as the total atmospheric error, as measured by 
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Figure 13. Strehl ratio at 8542 A for the data plotted in Figure fl2l These data have been 
smoothed by convolving them with a Gaussian kernel, with a standard deviation of 15. Al¬ 
though the Strehl ratio is at times higher than 0.2, when only the TTM is active, it is highly 
variable. 


KAOS. At around frame 3000, the system begins correcting for tip-tilt errors 
only. The residual error becomes the same as the total error minus the tip-tilt 
errors. Finally, at around frame 7000, full AO correction is applied, the residual 
error drops far below the total error lines. The seeing was quite good with, ro of 
around 17cm. 

Figure [13] Shows the corresponding Strehl ratios, calculated during the same 
sequence as that in Figure [T3] Although the seeing is good, the improvement 
allowed by correcting for tip-tilt errors alone is highly variable and much lower 
than that given by making the full correction. 
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5. Conclusions 


The performance of the new off-limb solar AO system has been demonstrated. 
This system has been shown to perform well, even under less-than-ideal seeing. 
It performs exceedingly well during good seeing. This is the first time to the 
knowledge of the authors that any group has used the light of solar prominences 
directly to measure the wavefront in an AO system. This allows for superior 
correction of solar prominence images, in the near infrared, when compared to 
what was previously available. 

The final Strehl ratio obtained by this system has been shown to exceed the 
predictions which was made. The Strehl ratio is high enough, under most circum¬ 
stances, to allow for the reconstruction of fully diffraction-limited images, at 8542 
A. This allows for spectroscopy and spectro-polarimetry at the diffraction-limit 
of the telescope. 

With the significant improvements allowed by the off-limb solar AO system, 
solar physicists will be able to further probe the mysteries of solar prominences, 
by deducing their properties at smaller spatial scales than previously possible. 
The authors believe that this will greatly improve humankind’s understanding 
of these phenomena and accompanying processes, such as CMEs. 
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